IPPP/09/03, DCPT/09/06, QMUL-PH-09-01 



Two-Loop Polygon Wilson Loops in J\f = 4 SYM 



O C. Anastasiou a , A. Brandhuber 6 , P. Heslop 6 , V. V. Khoze c , B. Spence b and G. Travaglinitl. 



O 

(N 

G 

1—5 

oo 

(N 



Oh 

> 

in 

(N 
(N 

(N 
O 
On 
O 



X 



a Institute for Theoretical Physics, ETH Zurich 
8093 Zurich, Switzerland 



b 



Centre for Research in String Theory 

Department of Physics, Queen Mary, University of London 
Mile End Road, London, El 4NS, United Kingdom 

Institute for Particle Physics Phenomenology, Department of Physics, 
Durham University, Durham, DH1 3LE, United Kingdom 



Abstract 



We compute for the first time the two-loop corrections to arbitrary n-gon lightlike 
Wilson loops in J\f = 4 supersymmetric Yang-Mills theory, using efficient numerical 
methods. The calculation is motivated by the remarkable agreement between the 
finite part of planar six-point MHV amplitudes and hexagon Wilson loops which has 
been observed at two loops. At n = 6 we confirm that the ABDK/BDS ansatz must 
be corrected by adding a remainder function, which depends only on conformally 
invariant ratios of kinematic variables. We numerically compute remainder functions 
for n = 7, 8 and verify dual conformal invariance. Furthermore, we study simple 
and multiple collinear limits of the Wilson loop remainder functions and demonstrate 
that they have precisely the form required by the collinear factorisation of the cor- 
responding two-loop n-point amplitudes. The number of distinct diagram topologies 
contributing to the n-gon Wilson loops does not increase with n, and there is a fixed 
number of "master integrals", which we have computed. Thus we have essentially 
computed general polygon Wilson loops, and if the correspondence with amplitudes 
continues to hold, all planar n-point two-loop MHV amplitudes in the M = 4 theory. 
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1 Introduction 



A surprising feature encountered in the study of supersymmetric gauge theories is 
the existence of an intriguing iterative structure in the higher-loop expansion of the 
Maximally Helicity Violating (MHV) scattering amplitudes in planar M = 4 su- 
per Yang-Mills (SYM) theory. This iterative structure was first discovered by Bern, 
Dixon, Kosower and one of the present authors (ABDK) studying collinear limits 
of maximally supersymmetric gauge theory amplitudes, and in the planar four-point 
MHV amplitude at two loops pQ. In the same paper, it was also conjectured that the 
same iterative structure should hold for two-loop MHV amplitudes with an arbitrary 
number of external legs. 

In a subsequent important development, Bern, Dixon and Smirnov (BDS) pro- 
posed an all-loop resummed formula for the n-point MHV amplitude, which they 
were able to confirm in an impressive three-loop calculation of the four-point ampli- 
tude pj. According to this conjecture, multi-loop amplitudes can be re-expressed in 
terms of the one-loop amplitude and four kinematic-independent functions of the 't 
Hooft coupling. One of these functions is the cusp anomalous dimension [3-5J, for 
which an all-order expression has been proposed in [6]. 

The ABDK/BDS conjecture was further investigated in several papers. In par- 
ticular, it was confirmed in a two-loop calculation for the five-point amplitude in [7] 
- a result which is particularly non-trivial since it implies a cancellation of certain 
parity-odd terms in the two-loop term of the logarithm of the amplitude^] It was also 
pointed out in [9] that the amplitudes of the /3-deformed M = 4 theory with real /3 
are identical to those of the undeformed theory (modulo an irrelevant, overall phase), 
and as such they will satisfy the ABDK/BDS iterative structure if the correspond- 
ing undeformed amplitudes do. Explicit expressions of the four-point amplitudes at 
four and five loops were also derived in [TU] and [UJ, respectively, and, for the four- 
dimensional cut-constructible part of the five-point amplitude at three loops in [12]; 
these expressions will allow for further tests of the BDS ansatz at four and five loops 
once the relevant integral functions have been evaluated to the necessary degree of 
accuracy in e. 

One of the key aspects of the ABDK/BDS conjecture is the appearance of the 
exponentiation of the one-loop result in the complete perturbative answer. In a 
remarkable paper [13], Alday and Maldacena succeeded in using the AdS/CFT cor- 
respondence to provide a string theory formalism to address scattering amplitudes at 
strong coupling. In particular, their calculation of the four-point amplitude repro- 
duced the strong-coupling limit of the BDS ansatz. It also provided a string theory 

lr The iteration for the parity-even terms had been proved earlier in [5]. 
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explanation for why planar scattering amplitudes at strong coupling exponentiate, 
through a semiclassical calculation. It was argued in [13] that the same exponentia- 
tion of [13] should hold not only for MHV amplitudes but also for non-MHV ampli- 
tudes, since the helicity dependence of the amplitudes in the prefactor is unlikely to 
modify the semiclassical exponent in the path integral. However, for the non-MHV 
case the exponentiation can only occur at strong coupling, and is not apparent in 
perturbation theory. 

The result of [13] suggested that the vacuum expectation value of a polygonal 
n-edged Wilson loop, evaluated this time at weak coupling, could be related to the 
perturbative n-point MHV amplitude in M = 4 SYM [T51HB]. This was confirmed 
in a one-loop calculation for n = 4 in [15] and subsequently for arbitrary n in [T6] . 
Drummond, Henn, Korchemsky and Sokatchev (DHKS) were later able to confirm 
this conjecture in a remarkable analytic calculation of the two-loop four-edged Wilson 
loop [T?j , followed by a semi-analytical calculation of the five-edged Wilson loop [32] • 

It was later argued [19] that the BDS ansatz may be incomplete, specifically for 
n-point amplitudes with n > 6 [20J. The authors of [21] have carried out an explicit 
calculation which shows that the BDS ansatz indeed needs to be modified in order to 
reproduce the two-loop term of the logarithm of the six-point amplitude. In a parallel 
development, the corresponding six-point lightlike Wilson loop was computed at two 
loops in [22j[23], and compared in [21j[23] to the parity-even part of the amplitude 
evaluated in [21] H The result of this analysis is that the MHV amplitude, stripped 
of the tree-level prefactor, and the Wilson loop are in perfect agreement (up to an 
additive constant) for the two-loop, six-point casejf] but there is an additional contri- 
bution compared to what the BDS ansatz predicts. This extra term, which we will 
refer to as the remainder function, will be one of the main characters of our paper. 

The possibility of having a nonvanishing remainder function was neatly explained 
in [TS] in terms of the anomalous conformal symmetry of the lightlike Wilson loop. 
In that paper the associated anomalous conformal Ward identities were derived, and 
it was shown that the BDS expression provided a particular solution to these Ward 
identities. At the same time, the anomalous Ward identities cannot uniquely deter- 
mine terms that are invariant under the conformal symmetry, and this leaves room 
for a conformally invariant remainder function [18]. For n < 5, the lightlike con- 
straints on the particle momenta restrict such conformally invariant contributions to 
just (kinematic-independent) constants. However, starting from n = 6 edges one can 
build functions of the conformally invariant ratios which are left undetermined by 
the Ward identities, and need no longer vanish. DHKS made the prediction therefore 
that, if the duality with Wilson loops holds, the remainder function should depend 

2 The full six-point amplitude at two loops has been presented in |24) . 

3 More accurately, there is a difference in the coefficients of the subleading f/e pole for the Wilson 
loop and the amplitude. We will come back to this point in Section 13721 
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on the kinematics only through cross-ratios. 

The dual conformal symmetry of the Wilson loop was also instrumental in sug- 
gesting that the S'-matrix of the Af = 4 theory should possess a dual superconformal 
symmetry [2t)J[2n], w hich is expected to be exact at tree level, and violated by an 
anomaly at the loop level. Indeed, it was later proved in [27] using a supersymmetric 
version [27H29] of the BCF recursion relation [30, 31] that the tree-level S'-matrix of 
the planar M = 4 theory is covariant under the dual superconformal symmetry. A 
solution of the supersymmetric recursion relation of [271429] was also presented in [32J . 
We also mention that the dual superconformal charges resurface as part of an infinite 
tower of charges coming from integrability of the dual AdS sigma model [331 154]. 

In this paper, we present the results of our study of the Wilson loop remainder 
function for arbitrary n. One important observation is that the structure of infrared 
and other integrable singularities of the diagrams which enter the Wilson loop calcu- 
lation does not change for n > 7. Therefore, with the same numerical routines we 
can evaluate Wilson loops for arbitrary n. We should note that our calculations were 
performed for Euclidean kinematics, but a generalisation to Minkowskian kinematics 
is possible. 

There are several interesting properties of the remainder functions which we have 
analysed. The first one is its conjectured dependence on the kinematics only through 
cross-ratios. In our study we have collected ample numerical evidence that confirms 
this expectation for n = 6, 7 and 8 points. The second aspect is the study of simple 
and multiple collinear limits of this function. Specifically, using universal factorisa- 
tion theorems [3"oT[3"8"] for scattering amplitudes one can predict [2T] the behaviour 
of the n-point amplitude remainder functions under simple collinear limits, namely 
H n — > 7^n-i • In this paper we will show that the appropriately defined Wilson loop 
remainder function 7?.Jf L has exactly the same collinear behaviour, namely 



Notice that no additional constant term appears on the right hand side of (11.11) . We 
have checked this numerically for n = 6, 7 and 8 sided polygon Wilson loops. Finally, 
one of the most important goals for the future is to find analytic expressions for 
the remainder functions. As a first step we have initiated a detailed map of these 
functions, in particular for n = 6, for a wide range of values of the cross-ratios. We 
were able to make intriguing observations for special values of the cross-ratios and 
lower dimensional slices of the kinematic parameter spaces. On general grounds, we 
expect the remainder functions to be transcendentality four functions of the conformal 
cross-ratios. However, even if we restrict the remainder functions to a one-dimensional 
slice of the parameter space, the space of transcendentality four functions is rather 
large and, hence, numerical methods are not sufficient to determine the remainder 
function. Clearly, new theoretical ideas, possibly from the AdS / CFT correspondence 
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or integrability, are needed in order to make progress in this direction. 

The rest of the paper is organised as follows. In Section 2 we review salient features 
of planar gluon scattering amplitudes in M = 4 SYM, their recursive properties at 
loop level, and the BDS all-loop ansatz. Furthermore, we introduce the amplitude 
remainder function, which is the difference between the full amplitude and the BDS 
ansatz, and is expected to depend only on the dual conformal cross-ratios of kinematic 
invariants. In Section 3 we set up the corresponding polygon Wilson loop calculations 
at two loops, and give a natural definition of the Wilson loop remainder function, 
which behaves under collinear limits in the same way as the amplitude remainder 
function. In Section 4 we present details about the diagrams and the corresponding 
Feynman integrals entering the calculation of arbitrary, lightlike, n-gon Wilson loops. 
In Section 5 we discuss the numerical evaluation of these Feynman integrals. In 
Section 6 we present a detailed, numerical analysis of the six-point remainder function 
including tests of dual conformal invariance, the explicit values of the remainder 
function at specific values of the cross-ratios, and various illustrative plots. In Section 
7 we start off with a discussion of the seven-point remainder function, and give explicit 
numerical results to illustrate our checks of dual conformal invariance and invariance 
under cyclic permutations and reflections of the external momenta. We then move 
on to discuss simple and multiple collinear limits of two-loop amplitudes and Wilson 
loops. Specifically, we present evidence that in the simple collinear limits the seven- 
point remainder function becomes equal to the six-point remainder function. Finally, 
in Section 8 we present a similar analysis for the eight-point Wilson loop and briefly 
discuss the generalisation to arbitrary n. 



2 Planar amplitudes in J\f = 4 super Yang-Mills 
and the ABDK/BDS ansatz 

The infinite sequence of n-point planar MHV amplitudes in TV = 4 SYM has a re- 
markably simple form. At any loop order L, the amplitude can be expressed as the 
tree-level amplitude, times a scalar, helicity-blind function 

=AT e Mi L \ (2.1) 

At one loop, the function Ain^ is simply a sum of two-mass easy box functions 
F 2me [39], with coefficient equal to one: 

M (i) = J2F 2 ™ e (p,q,P,Q) . (2.2) 
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where 
and 



In [T], ABDK discovered an intriguing iterative structure in the two-loop expansion 
of the MHV amplitudes at four points. This relation can be written as 

Mf{e) - \{M^\e)) 2 = f^(e)M?(2e) + C« + 0(e) , (2.3) 

f {2 \e) = -C 2 -C 3 e-C4e 2 , (2.4) 

= -I C | . (2.5) 

In pp, it was conjectured that f l2.3p -( l2~5|) should hold for two-loop amplitudes with 
an arbitrary number of legs - a conjecture which was consistent with an explicit 
evaluation of the universal two-loop splitting amplitude. 

Building upon the iterative relation of pQ, and the known universal infrared be- 
haviour of gauge theory amplitudes |4T)144"9"] . BDS proposed a resummed, exponen- 
tiated expression for the scalar function Ai n . In the same paper, this conjecture 
was checked in a three-loop calculation in the four-point case. Specifically, the BDS 
conjecture is expressed as [2] 

oo oo 

M n := l+J2* L Mi L \e) = exp[Y,a L (f L \e)M£\Le)+C^+Ei L \e))] , (2.6) 

L=l L=l 

where a is the loop-counting parameter. In the conventions of [2], this is defined as 
a = [# 2 iV/(87r 2 )](47re- 7 ) e . Here f (L \e) is a set of functions, 

f {L \e) := ti L) + f[ L) e + ti L) e\ (2.7) 

one at each loop order, which appear in the exponentiated all-loop expression for 
the infrared divergences in generic amplitudes in dimensional regularisation [17] (and 
generalise the function in ( 12. 3p ). In particular, = 7^/4, where jk is the 
cusp anomalous dimension, 

oo 

lK {a) =J2 aL > Ik = 4 , Ik =-4(2, (2.8) 

L=l 

related to the anomalous dimension of twist-two operators at large spin [H]. The 
appearance of the cusp anomalous dimension characterises the relation between the 
infrared divergences of scattering amplitudes and ultraviolet divergences of Wilson 
loops with cusps which was originally discussed in QCD in [51 14*3] . 

The 0(e) term in (12. 7p is related to the so-called collinear anomalous dimension 
G, f[ L) = (L/2)GW, 

oo 

G(a)= Y,* L G {L) , G^ = -( 3 , (2.9) 



and /j = —(4- In particular, can already be found from simple collinear limits 
of the amplitudes. For future reference we also define 

00 

C{a) = flL ° {L) • ( 2 - 10 ) 

Importantly, the constants C^ L \ /q , fx and 011 the right hand side of ( 12. 6 p do 
not depend either on kinematics or on the number of particles n. On the other hand, 
the non-iterating contributions depend explicitly on n, but vanish as e — > 0. 

BDS also suggested a resummed expression for the appropriately defined finite 
part of the n-point MHV amplitude, 

Tn = e^ DS , (2.11) 

where 

F n BDS (a)= ^(a)F«(0) + C(a). (2.12) 

The quantities Jk{o) an d C(a), are given in (I2.8P and (12. 10 j) : the entire dependence 
on kinematics of the BDS ansatz enters through the finite part of the one-loop box 
function, Fn(0). Explicitly, one has [2] 

-Mi 1) ( e ) = -^E(-^j +^\e), (2.13) 
1 " 

1=1 

where 

[n/2]-l / [r] \ / [r] \ 

9n,i = - £ In I -^J hl ( I#T J + °^ + L - + 2 C2 ' (2 ' 15) 

and tf := (p» + • • • +Pj+ r _i) 2 are the kinematical invariants. The explicit forms of 
the functions D n i and L n ^ depend on whether n is odd or even. For n = 2m one has 

m-2 / t [r k [r+2] \ 1 / t im-l] t [m+l]\ 

D 2 m,i = ~ g Li 2 f 1 - ^^xff+i] ) ~ 2 Lla 1 1 " ' t M t M J ' (2 ' 16) 



L 2mi = _ii n i* h 4t 

\ L i+m+l / \ L i+m, 



1 



whilst, when n = 2m + 1, 



m-l / t [r] t [r+2] 
r=2 V h 



J 2m+l,i 



1 / j.L"'J \ / -fl m l 

In/ 



2 \ _+H / I +M 

' i+m+l / V i+m , 



The case n = 4 is special; in this case the finite remainder is given by 

F 4 (1) (0) = ^ln 2 ^) +4C 2 , (2.18) 



2 \t 

which is the finite part of the zero-mass box function plus a constant shift 



2.1 The amplitude remainder function and n-point cross- 
ratios 

Beyond five points, and starting from two loops, the ABDK/BDS ansatz (12. 3p needs 
to be modified by the addition of a remainder function lZ r , 



M^(e) - \(M%\e)) 2 = f^(e)M^(2e) + + <R n + (9(e) . (2.19) 
We now move on to characterise it. 

To begin with, we recall one of the important properties of the ansatz, namely 
that it already incorporates the correct simple collinear limits of the amplitude for all 
n [HE]- We will review this in Section I7T2"| but we would like to anticipate here one 
important consequence of this, namely the fact that the remainder function must have 
trivial collinear limits. With the definition given above of the remainder function, one 
expects that under a simple collinear limit [21] 



n n -> ft n _i , (2.20) 
where no constant term can appear on the right hand side of (I2.20p . 

An important advance was made in [18J, where it was realised that the BDS ansatz 
is a solution to the anomalous Ward identity for the Wilson loop associated to the 
dual conformal symmetry. As is by now common, one introduces dual (or region) 
momenta x±, ■ ■ ■ ,x n [T3l[50] and defines the particles' momenta as 

Pi := Xi - x i+1 , (2.21) 

(with the identification x\ = x n+ i), which satisfy n on-shell relations (x^ — Xj + i) 2 = 0, 
i = 1, . . . , n. The dual conformal group then acts on the dual momenta pTJlIE]. It is 
important to notice that acting with conformal transformations on the dual momenta 
x does not endanger momentum conservation, which is automatically satisfied once 
the momenta are written in the form ( 12.211) . 

Obviously, adding to the ABDK/BDS ansatz any arbitrary function of the con- 
formally invariant cross-ratios 

, (2.22) 
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preserves conformal invariance and hence would give another solution to the same 
conformal Ward identity. With this in mind, it was argued in [18] that the remainder 
function, if non-vanishing, should depend on the kinematics of the scattering only 
through cross-ratios. It is therefore important to examine how many independent 
cross-ratios one can build at n points. 

Starting from a set of n arbitrary points, one can define n{n — 3)/2 cross-ratios. 
This number is reduced to n(n — 5)/2 if one imposes the n on-shell conditions for 
lightlike momenta - notice that this is the same as the number of two-mass easy 
box functions which could potentially appear in a colour-ordered Yang-Mills n-point 
amplitude. If one considers four-dimensional external momenta we have additional 
constraints that the Gram determinant of any five of them should vanish. With the 
Gram determinant constraints, the number of possible on-shell cross-ratios is reduced 
to [18] 3n — 15 for n > 5. Incidentally, we note that the number of on-shell cross-ratios 
one can construct from n points in a .D- dimensional space with D > n — 1 is also 
equal to n(n — 5)/2 [TB]. Discarding the Gram determinant constraints is therefore 
equivalent to considering the external momenta as defined in a .D-dimensional space 
with D > n — 1. 



Taking this apparent coincidence between the number of cross-ratios and of two- 
mass easy boxes more seriously, we define the independent cross-ratios we will use to 
parametrise any n-point remainder function as 

«« == 3^ ■ (2-23) 

In [16] it was shown explicitly how a finite one-loop Wilson loop diagram, where a 
gluon propagator connects the momenta Pi = Xi — x i+ i and pj = Xj — Xj + i, reproduces 



the finite part of the two- mass easy box function with kinematic invariants s = x~j, 
t = x? +1 j +1 , P 2 = xfj +1 , Q 2 = x 2 +l y The choice of kinematic invariants appearing in 
( I2.23P precisely matches those appearing in the corresponding box, and tty is of the 
form P 2 Q 2 /(st) (see Figure H]fl 

Two comments are in order. Firstly, we observe that the ratios in (I2.23P are the 
same as those entering the functions D and L defined in (|2.16p . (I2.17p . which appear 
in the BDS ansatz (I2.12p . Secondly, we mention that a generic cross-ratio, i.e. one 
of the form x 2 jX 2 m / (a^a;| m ) for generic k, I , can be written as a product of cross- 
ratios of the form (I2.23P for arbitrary n. If we assume i < I and j < m this is simply 
given by 

2 2 '—1 rn—1 

= n n «~ ■ p- 24 > 

ij Im r=i s=j 



1 Note also that a basis of off-shell cross-ratios can be given similarly, simply by allowing also the 
one-mass box functions as well as the two-mass easy box functions. 
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Figure 1: On the left we represent the one-loop Wilson loop diagram which gives the 
finite part of the two-mass easy box function with massless momenta Pi,Pj [16J. On 
the right we represent the corresponding cross-ratio Uij, the red dashed lines depicting 
the factors xf +1 j +1 , xf +1 j, xjj +1 in the definition ofuy in ( I2.23p . 

The two different countings of cross-ratios mentioned above - with and without 
Gram determinant conditions taken into account - predict that no cross-ratios can be 
written for n = 4,5, whereas at six points they allow for three independent harmonic 
ratios. These could be chosen to be 

2 2 2 2 2 2 

rp*" nn^ fY* rp" rp" rp" 

_ x 31 x 46 ._ _ X 15 X 24 ._ _ x 26 x 35 ._ ( 9 r\ 

x 36 x 41 • i 14 x 25 x 25 x 36 



In [21], [23] it was verified for several kinematical configurations that the six-point 
remainder function indeed depends on the kinematics only through the three cross- 
ratios in (I2.25p . as predicted by dual conformal invariance. Furthermore, the six-point 
remainder function is a symmetric function of the three cross-ratios (I2.25P [22] . It 
can easily be shown that any permutations of the three arguments of the remainder 
function corresponds to a cyclic relabeling of the dual momenta plus possibly a rever- 
sal of their ordering, which clearly leaves the Wilson loop unchanged. This property 
was checked numerically in 



Interestingly, in [23] it was also shown numerically in a few examples that two 
kinematical configurations which have the same cross-ratios, but differ in that one 
respects the Gram determinant constraint and one violates it, give rise to the same 
numerical values for the remainder function. In this paper we perform explicit calcu- 
lations of lightlike Wilson loops at six, seven and eight points. Starting from seven 
points, we find that there is a different number of cross-ratios depending on whether 
one implements the Gram determinant constraint or not. 

In practice, in the following we will discard the Gram constraint altogether, and 
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work with unconstrained kinematics. This turns out to be a particularly efficient way 
to generate kinematical points - including cases where n is odd. At seven points, 
we will therefore consider seven (rather than six) cross-ratios of the form (I2.23p . At 
eight points, we will consider twelve cross-ratios (rather than nine). In all cases, we 
have performed extensive numerical checks proving that the remainder function only 
depends on kinematics through the expected cross-ratios. 



3 Wilson loops and scattering amplitudes 



The Wilson loop we consider in this paper is purely bosonic, and its expression is 
given by 



W[C n \ := TrPexp 



ig (f) dr i"(r)A M (x(r)) 



(3.1) 



The particular closed contour C n we consider is the lightlike n-edged polygonal contour 
introduced in [13]. It is obtained by attaching the momenta of the scattered particles 
Pi, . . . , p n one after the other, following the order of the colour generators in the colour- 
ordered scattering amplitude. The resulting contour is closed as YH=iVi — 0, and 
the positions of the vertices are given by the dual momenta coordinates, introduced 
earlier in ( I2.2ip . 

Calculations of (13.1 p at one loop were performed in [15] and [16] , where agree- 
ment was found with the expression of the scalar function in (12. 2p appearing in the 
corresponding one-loop MHV amplitude. In the following we will discuss the basic 
ingredients needed to perform a two-loop perturbative calculation of the Wilson loop. 



3.1 Perturbation theory setup 



The calculation of the Wilson loop at higher loops is simplified if one makes use of 
the non-abelian exponentiation theorem - This theorem allows one to write 

the result of the vacuum expectation value of the Wilson loop as an exponential, and 
gives a practical rule to calculate the exponent. We represent the Wilson loop as 

oo oo 

(W[C n \) : = 1 + $>'W« := exp , (3.2) 
i=i i=i 
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and, in this paper, our main focus is on the evaluation of the two- loop term in 
the exponent. In terms of the Wilson loop coefficients this is obtained as^| 

^ = Wf ) - \ ( W Mf . (3.3) 

The non-abelian exponentiation theorem has been used widely in several Wilson loop 
calculations, see for example [53j[5l], and more recently [T71EEE1I23]. To briefly il- 
lustrate its application, we first consider the calculation of a Wilson loop vacuum 
expectation value in an abelian theory. In this case, it is not difficult to see that the 
perturbative series reorganises itself into the exponential of the one-loop correction, 
i.e. the corresponding abelian result is given by a formula like (13. 2 p with Wq ED = 0, 
for any I > 1. In the non-abelian case, parts of the result of the diagrams contribute 
to the exponentiation of the one-loop result, but there are additional contributions 
which correct order by order in perturbation theory the one-loop term in the expo- 
nent. In brief, the rule for calculating the complete exponent [511122] * s t° restrict to 
those parts of the diagrams which give a "maximal non-abelian colour factor". At 
two loops, this turns out to be equal to CfCa where Cp := C^ir) is the Casimir 
in the representation r of the Wilson loop, and Ca '■— C*2(G) is the adjoint Casimir |f] 
For SU(N), one has Cf = {N 2 — l)/(2N) for the fundamental representation, and 
C A = N. 

As a simple example, consider the two-loop contribution to a cusp diagram arising 
from diagrams containing only propagators. The contribution from a ladder diagram 
produces the colour factor r Jr(T a T a T 6 T 6 ) = dpCp, whereas the cross propagator 
diagram, represented on the left hand side of Figure [7J contains the colour factor 
Tr(T a T 6 T a T b ) = cIfCf{Cf — 1/2 Ca)- According to the non-abelian exponentiation 
theorem, we only have to consider the term —(1/2) CfCa from the cross propagator 
diagram, and discard the remaining diagram altogether (which has already been taken 
contributed to the exponentiation of the one- loop correction). 

As a final remark, we would like to observe that the diagrams needed to calculate 
these maximally non-abelian corrections are simpler (and fewer) than those needed for 
the full Wilson loop, however the technical difficulties in obtaining the final integrals 
in analytic form are typically comparable. 

We now move on to describing the basic ingredients of any Wilson loop perturba- 
tive calculation. The first one is the gluon propagator which, in the Feynman gauge, 

2 Our Wilson loop conventions are summarised and compared to those of [25] in Appendix A. 
3 We notice that, in order to be properly normalised, the Wilson loop in (|3.ip should be divided 
by the dimension of the representation dp := d(r). 
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is given by A M! ,(x) := r] lxv A{x), where 



7r"a (D \ 1 
47T V2 /(_ x 2 +i£ w-i 



vr £uv iYl-e 



47r 2 (-x 2 + ieY- evv 



where D = 4 — 2e uv . The Wilson loop is gauge invariant, therefore we can pick any 
gauge we like to compute its expectation valued 

At two loops, we will also have Feynman diagrams where the gluon three-point 
vertex contributes. The basic structure to know is therefore the Wick contraction of 
three gauge fields with a three-point vertex, 



Wick 



A" 1 ' ax {xi)A" 2 ' a2 {x 2 )A'" 3 ' a3 {x 3 ) Jd D zTr(d tl A 1/ [A' x ,A u }yz)j = -iC F f 



eaia 2 a 3 



where 



G(xi, x 2 , £3) := Jd D z A(xi — z)A(x 2 — z)A(x 3 — z) (3.6) 



64vr D -1 7=\ 1 it 1 i a i a ^i2 + ai«3^?3 + «2a3^i 3 )^ 3 ' 

where we have used (13.41) and xjj = (xj — Xj) 2 . The evaluation of the right hand side 
of (13 M in various cases, specifically when x\ = x\ = x\\ = 0, has been carried out 
in 



Finally, we notice that the colour factor associated with gluon three-point ver- 
tex diagrams, obtained after contracting with a trace of three colour generators, is 

^rpaxrpa^rpaa^faia^as = (\/2)d F C F C A - 



3.2 The Wilson loop remainder function 

We define the n-sided Wilson loop remainder function TZ^ h in complete analogy with 
the amplitude remainder function introduced in (12.191) . as 

= f%{(e) W M(2e) + C« + 72™ , (3.7) 

4 The advantage of considering different gauges, still belonging to the class of Feynman-'t Hooft 
gauges, has been discussed recently in [55]. A different possibility would be to pick the lightcone 
gauge. This gauge has been used for Wilson loop calculations in [56] . 
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where e = — e uv . We have added a subscript WL to distinguish quantities relevant 
for the Wilson loop from the corresponding amplitude expressions. In particular, 

/wl( e ) := /o (2) + /fwL 6 + /iiwi/ 2 > ( 3 - 8 ) 

where is the same as in (ED, whilst /fi L = G%1 = 7C 3 [53], and the third term 

/ 2 WL is a so far undetermined constant. Similarly, the constant C^ L in (I2.19P has 
also not been fixed yet. We will shortly determine these two constants. 

A few comments are in order. 

(2) 

Firstly, we notice the already mentioned discrepancy between the coefficient G ei ^ 
of the subleading 1/e pole in the Wilson loop and the corresponding coefficient 
appearing in (12. 9p on the amplitude side. This discrepancy has been examined and ex- 
plained in [FT]. We note that this discrepancy cannot be reabsorbed into a (kinematic- 
independent) redefinition of the Wilson loop renormalisation scale /iwl alone |^ 

Secondly, we would like to point out that if the correct determination of the 

(2) (2) 

constants / 2 wl anc ^ C*wl * s implemented also for Wilson loops, we expect the Wilson 
loop remainder function to have precisely the same collinear limit as its amplitude 
counterpart, namely 

K^^K™, (3.9) 
with no extra constant term on the right hand side of (13. 9p . 

In order to determine C^ L and f% % L , and later be able to check ( 13. 9p . we proceed 
as follows. Firstly, we recall that conformal invariance guarantees that the four- and 
five- point Wilson loops satisfy an ABDK/BDS-like ansatz [18J. This implies that 
the remainders TZ^ L and TZ^ L cannot depend on kinematics and must be constant. 
On the amplitude side, these remainder functions are known to vanish. Thus we also 
choose 

^2'\ (2) _____ 

in (13. 7p and then determine C WL an d / 2 wl fr° m solving (13. 7p for n = 4, 5. 

Notice that in writing the Wilson loop ABDK/BDS ansatz, it is crucial to use 
the one-loop Wilson loop, and not the one-loop amplitude. The two are equal to all 
orders in e up to their normalisation [TB] . More concretely, for the amplitude we have 
bddk where Ai^ n bddk * s the one-loop amplitude in the normalisations 

of 135] , and 



gra + ^q-.) 

2 r(l-2e) V ; 



5 See also the discussion in [33] (version 3). 



13 



L5DDK- On * ne other hand, 



This leads to M ( n } = e er T(l + e)(r(l - e) 2 /r(l - 2e))M$ 

for the Wilson loop, we hav^ w$p = e e T(l + e) ■M^bd'oK' This leads to the following 
correspondence between the Wilson loop and the amplitude at one loop, 



r(l-2e) (1) 
r 2 (i — e) 



;i + C 2 6 2 )M 1 ) + 0(e 



- + 0{e) . (3.12) 



At one loop, the four- [15] and five-edged Wilson loops [16] are thus given by 



(i) 



w 

2 

8=1 

and at two loops 




(3.13) 
(3.14) 



(2) 



(2) 



5 

E 

i=l 



J2]' 



/// 



-2c 



-2e 



7T~ 



48e 2 




7T 

72 



(3.16) 



We note that in ( 13.15)) we have used the results of our two-loop calculation of the four- 
point Wilson loop to correct the constant term in the corresponding result of [T7] 

We can now uniquely rewrite (13. 15j) and ( I3.16P in an ABDK/BDS form as 

wf(e) = f^(e)w?(2e) + C^, (3.17) 

(3.18) 

(3.19) 
(3.20) 



w 



(2), 



e) = f$ L (e)w^>(2e) + a 



.(2) 
'WL ' 



where 
and 



(2) / 
WlA 



-Ca + 7(3 e - 5C 4 e^ 



a 



(2) 
WL 



-c 2 



The 0(1) and 0(e) coefficients of f^ L (e) had already been determined in [TT] . Inter- 
estingly, the constant C^l turns out to be the same as the constant in (12. 5p for 
the amplitude. 



6 In the following formulae we employ the redefinition of the renormalisation scale in (|4.3|) . 
7 This discrepancy has also been noted independently by Marcus Spradlin, whom we thank for 
discussions on this point. 
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Finally, let us now compare the definition of the remainder function 7cg given in 
(13. 7p with that of DHKS [22i[23] (see also Appendix A). First, we write the two-loop 
term in the form 

n / 2 \ -2e (2) (2) 
w (2) = £ (_^\ ^ + ^ + F (2) + 0(£) j (3 21) 

(2) 

where Fn is finite as e — > 0. 



From ( I3.15P and ( I3.16p . one has 

"3 = £ , = , (3.22) 

and comparing (I3.2ip with (13. 7p . using (I3.12p and ( I2.13p . we obtain 

if } = \lP FM(0) + C%{ + n^+n n -. (3.23) 



The DHKS finite remainder function is then defined as 

7^ HKS : = - F n BDS(2) , (3.24) 

where Fn DS ^ is the two- loop contribution to F^ BS in ( 12.12[) . 

if DS(2) = \iPf£\0) + CV. (3.25) 

Thus we find that our finite remainder defined in (13.71) and the DHKS definition (13.241) 
are related by a constant shift, 

n^ = K^-n^-. (3.26) 

We have checked that our Wilson loop remainder function 1Z^ L satisfies ( 13. 9p under 
a collinear limit. For n = 6 this amounts to 1Z^ L — > lZ^ h = 0. This would imply 
that 

ftDHKS _^ n 4/ g = 12 .i76l... , (3.27) 

which is completely consistent with the results of [22], where this was computed 
numerically as cw := 12.1756 with accuracy of the order 10~ 3 . In the following 
sections we will show that (13.91) also holds for n = 7, 8. 



We can now express the statement of the duality between Wilson loops and am- 
plitudes as an equality of the corresponding remainder functions defined in (13.71) and 

dsn, 

TZ n = K™ L . (3.28) 
Notice that no additional constant term is allowed on the right hand side of (I3.28p . 
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4 Summary of the diagrams entering the Wilson 
loop at two loops 



In this section we summarise the expressions for all the diagrams entering a generic 
two-loop Wilson loop calculation. There are five main ingredients to the two-loop 
Wilson loop calculation at any number of edges, n. We call them the "hard diagram" , 
the "curtain diagram", the "cross diagram", the "Y diagram", and the "factorised 
cross diagram" , see Figures [2Hni 

In the following we summarise the final expressions for the integrals corresponding 
to these diagrams; derivations are outlined in the Appendices. The entire two-loop 
contribution to the logarithm of the Wilson loop is assembled in terms of the indi- 
vidual building blocks in Section 14.61 

In all expressions of the diagrams we will write in the next sections, a factor of 

2 2 

C := 2a 2 /i 4£ [r(l + e)e 7e ] = 2a 2 // e (l + y e 2 ) + £>(e 3 ) , (4.1) 
will be pulled out, where we have defined the coupling 

and the scale, /i 2 , is given in terms of the Wilson loop scale as 

Mwl : = . (4.3) 

This factor will be reintroduced when the diagrams are reassembled into the complete 
Wilson loop (I4.12p in order to match the conventions of [23] and facilitate comparisons 
with their results. 



4.1 The hard diagram 

The hard diagram is depicted in Figure [2} and is given by the integral: 

fH(Pl,P2,P3;Ql,Q2,Q3) 

1 I! ( n *o jf ( n **) ^ - 1 <«) («^«3)--^ uv 



(4.4) 
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where the functions T> and M are given in ( 1B.2j) . ( IB. 51) . The momenta pi are massless 
pf = and the Qi can be massive. They are further constrained by momentum 
conservation 

Pi +P2+P3 + Q1 + Q2 + Q3 = . (4.5) 




Figure 2: The hard diagram. 



4.2 The curtain diagram 



The curtain diagram is represented in Figure El and its expression is given by the 
integral 

fc{vuP2,Pz\QiiQ2iQz) 



JO V „■_-, J Jl-n 



2 Jo \ =1 ' Ji-n [- 2(p 1 Q 3 )a 1 - 2(p 1 p 2 )a 1 r 2 - 2(p 2 Q 3 )r 2 - Qf] 1 

x — . (4.6) 

[ - 2(piQ 2 )n - 2( PiP3 )t 1 t 3 - 2(p 3 Q 2 )r 3 - Q|] 1 euv 

The and are constrained as in (14. 5p . 
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/'I 



Ml 




Figure 3: The curtain diagram. 



4.3 The cross diagram 



The cross diagram is represented in Figure HI and is given by 



fx(Pl,P2',Ql,Q2) 

1 per i r-T2 

:= — - / d(JidT 2 I dr\ J da 2 
^ Jo Jo Jo 



l-ei 



( - 2(p 1 p 2 )cria 2 - 2p 1 Q 2 cr 1 - 2p 2 Q 2 cr 2 - Qj) 

(P1P2) 

( - 2(p 1 p 2 )r 1 T 2 - 2p 1 Q 2 r 1 - 2p 2 Q 2 r 2 - Qf) 1 

(4.7) 

Again the Pi are massless and the Qi massive and momentum conservation is imposed, 

Pi+P2 + Qi + Q2 = . (4.8) 
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Pi 




l>2 



Qi 



Figure 4: The cross diagram. 
Pi Pi 




Pi 



Qi 



(h 



Figure 5: The Y diagram together with the self-energy diagram. 



4.4 The Y diagram + self-energy diagram 



The Y diagram, to which we also add (half of) the self-energy diagram J§ represented 
in Figure [51 is given by the following integral, 



fY{pi,P2]Qi,Q2) 



i ,i 

da / dr\dr2 

n JO 



Pi-Pa 1 r(l-2e uv ) 
8 e uv r 2 (l-e uv ) 

(J- £uv (l - cr) _euv 



(4.9) 



-Qi ~ 2(Q lP2 )r 2 - 2(Q lPl )ar 1 - 2(p 1 p 2 )ar 1 T 2 ) 1 ^v 
cr _euv (l - a) _euv 



-Ql - 2(Q 2 p 2 )r 2 - 2(Q 2 p 1 )aT 1 - 2(p 1 p 2 )ar 1 r 2 ) 1 - 



2c, 



3 The other half of the self-energy accompanies the "upside-down" Y diagram. 
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4.5 The factorised cross diagram 



Pk 



pi 



Figure 6: The factorised cross diagram. 
This is given as —1/2 times the product of two one-loop diagrams 

-i/2fp(Pi,Pj] Qji, Qij)fp{Pk,pi] Qik, Qki) ■ 



(4.10) 



The one-loop terms come from a diagram involving just a single propagator, and are 
simply the finite part of the two-mass easy box function [16] 



f P (j>,q;P,Q) = - 



il„g>(£) + Li 2 (l-^) + Li 2 (l-^) + Li 2 (l-^ 



+ Li, (1-^1 !• f I 



P 2 Q 2 



(4.11) 



where s = (P + p) 2 and t = (P + g) z 
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4.6 The complete Wilson loop at n points 



The logarithm of the complete n-sided Wilson loop is given by the sum over all 
diagrams, 

c \ S [fH(pi,Pj,Pk]Qjk,Qki,Qij) + fc(pi,Pj,Pk]Qjk,Qki,Qij) 

l<i<j<k<n 

+ fc(Pj iPkiPi] Qkii Qiji Qjk) 4" fc\PkiPiiPj] Qiji Qjki Qki) 

+ [fxipuPj] Qji, Qij) + fy(pi,Pj] Qji, Qij) + fr(pj,Pi] Qij, Qji) 

\<i<j<n 

+ Yl (- 1 / 2 )fp(Pi,Pj;QjhQij)fp(Pk,pr,Qik,Qki)>, (4.12) 

l<i<k<j<l<n ' 

where the first sum is over all sets of three non-equal legs k, the second sum is 
over all sets of two non-equal legs i, j, and the third sum over all sets of four non-equal 
legs. Here we have defined 

Qij = Pi+i + Pi+2 + ■■■ Pj-i , (4.13) 
and C is the factor we pulled out (see (14.1 1) ). 

Note that the singular properties of these integrals depends on whether Qi = or 
not (i.e. whether legs are adjacent or not). For example fn has a l/e uv 2 singularity 
if Qi = Q2 = 0, Q3 7^ 0, a l/e uv singularity if Q\ = 0, Q2)Qz 7^ 0, and is finite as 

c uv 

In the four-point case, for example, f)4.12p leads to 
wf ] =c|/ff(pi,p2,p3;0,p4,0) 

+ /c(Pl,P2,P4;P3,0,0) + fc(Pl,P2,P3-,0,Pi,0) + fc(PuP3,P4,0,0 } p 2 ) 
+ ^fx(Pl,P3-,P4,P2) + fx(pi,P2\Pl+P2,0) 

+ fY(pi,P3-,P4,P2) + /y(pi,P2;pi+P2,0) + Mp2,Pi;0,pi+p 2 ) 

+ (- 1 /8)fp(Pl,P4,P2)fp(P2,P4,Pl,P3) 

+ cyclic permutations of (pi,P2,P3,P4) | (4-14) 

(the factor of 1 /2 in front of fx and the extra factor of 1/4 in from of the factorised 
cross is to account for the double counting of diagrams when summing over cyclic 
permutations). Of course everything should only depend on s = (pi + p 2 ) 2 and 
t = (p! +p 4 ) 2 . 
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4.7 Cusp diagrams 



The formula for the exponent of the Wilson loop at two loops (14.121) includes all 
contributing diagrams. A subset of these diagrams involve only two consecutive edges 
and are known as cusp diagrams. These are given by 



n 

C ^2 (fx(Pi,Pi+i,Qi+ii,0) + /y(pi,Pi + i;Q i+ ii,0) + fY(Pi+i,Pi;0,Qi+ii)^) , (4.15) 
i=i 



and are shown in Figure [3 



The final result for the two-loop correction to the cusps is given by 



C f Ca 



,2 \ 2 



47T 5 



r(i + e)vH 



1=1 



X 



-2c 



ii+2 



fJ-WL ) 



/cusp(^) 



(4.16) 



where 



cusp V 



1 1 

28e^ 



r(l + 2e)r(l-e) 

r(i + e) 



7T 



^ IT + tStT + °( e ) ■ ( 4 - 17 ) 



48e 2 8e 160 



This way, one can rewrite (14 .16}) using the redefinition of fi (14. 3 1) as 



i=i 



X ii+2 



-2c 



r(l + e)e 



7C 



/cusp(^) 



(4.18) 



Since 



T(l + e)e 



1 2 



7C 



7T 



cusp V 



Cs 

48e 2 8e 



_____ 

720 



0(e) 



we obtain that the contribution from all cusps is therefore 



(4.19) 



2a 2 



,r 



ii+2 



-2c 



7T 



48e 2 8 



_____ 

120 



(4.20) 



Equation (I4.20p is exactly equal to the corresponding cusp results in [23] (after con- 
sidering that their e is the ultraviolet parameter). 



5 Evaluation of Wilson loop diagrams 



An intriguing property of the n-point polygon Wilson loop parameterisations from 
Section H] is that they are valid for an arbitrary n. The infrared properties of the 
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Figure 7: Maximally non-abelian Feynman diagrams contributing to the two-loop cusp 
corrections. The second diagram appears with its mirror image where two of the gluon 
legs of the three-point vertex are attached to the other edge; these two diagrams are 
equal. The blue bubble in the third diagram represents the gluon self-energy correction 
calculated in dimensional reduction. 

integrals needed to evaluate the diagrams may change from more to less divergent, 
however all cases can be obtained from the same starting expressions involving para- 
metric integrals with five dimensions at most. This is to be contrasted with the 
Feynman parameterisations of Feynman amplitudes, which require a different repre- 
sentation with additional parameters as the number of external legs n is increased. 

Our goal is to construct an algorithm which evaluates the diagrams of the previous 
section in complete generality, for an arbitrary number of points, n. This has been 
achieved for up to n = 6 in the literature [25]. Given the fact that the number 
of Feynman parameter integrations and the number of distinct "master" functions 
required for an evaluation at different values of n is independent of it, one could aim 
for solving this problem in complete generality. It is very exciting that Wilson loops 
and M = 4 SYM planar MHV amplitudes are very likely to be dual to each other. If 
this is proven correct, solving the problem of calculating two-loop n-gon Wilson loops 
also provides a solution to the problem of evaluating planar two-loop amplitudes with 
an arbitrary number of legs. 

A fully analytic evaluation of the master functions in arbitrary n-point Wilson 
loops appears to be a formidable task. A more viable and practical approach is to 
evaluate these integrals numerically. Many of the required integrals develop diver- 
gences when e = 0, with 1/e 2 poles at most. We use the programs developed in [58H6T] 
for evaluating generic Feynman diagrams in order to automatically cast these integrals 
as Laurent series in e. The coefficients of the series are multidimensional integrals, 
however they are free of singularities in the dimensional regularisation parameter, 
and we can evaluate them numerically using well established stochastic integration 
methods [62] . 

We notice that we have evaluated the finite "hard diagrams" (which appear for 
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the first time at n > 6 points), using an alternative approach; this is possible because 
the diagram is free of singularities in the limit e — > 0. We recall that a one-loop scalar 
triangle can be written using Feynman parameters as 

Tria(D; u u u 2 , u>) := J j4j^yn [{k + Pl) ^ [{k + Pl + p 2 )T ^ 

3 3 _1 V2 — 1 V3 — 1 



; r(^)r(v 2 )r(^ 3 ) 7 ^ % ia2P ? + + a 



[a\a 2 p\ + a 2 a 3 p 2 + a^aip^y 123 2 

with p 3 = — pi —p 2 and z/ 123 := V\ + z/ 2 + z/ 3 . Interestingly, this is the same function as 
the triple-gluon vertex of (I3.6p . In the absence of divergences, as occurs in the hard 
diagrams with nonvanishing Qi, Q 2 and Q3, we can set D = 4 (equivalently e = 0) in 
(13. 6p . By comparing (13. 6p and (15. ip . we observe that the triple-gluon vertex is a sum 
of one-loop scalar triangles in D = 6 dimensions, with powers of propagators taking 
the values Vi = 1,2. 

We exploit the mapping of the vertex of (13. 6 p onto familiar one-loop integrals to 
simplify it before we insert its expression into the hard diagram. Using an automated 
reduction program [53] , we express the finite part of the vertex as a linear combination 
of triangle Tria(6 - 2e; 1, 1, 1), and bubble Tria(4 - 2e; 0, 1, 1), Tria(4 - 2e; 1, 0, 1), 
Tria(4— 2e; 1, 1, 0) master integrals. Notice that we have re-introduced the dimensional 
regulator e because these master integrals are all divergent as e — > 0. However, we 
can set e = after we substitute the analytic expressions for the bubble master 
integrals, and a Feynman representation derived from (15. ip for the triangle master 
integral in six dimensions. An alternative basis of master integrals can be obtained 
using dimensional shift identities [64"||65]. where the triangle master integral is also 
in four dimensions. However, this choice is inconvenient for numerical evaluations, 
because the dimensional shift generates a Gram determinant in the denominator 

1 



a 3 (p?,p 2 ,p1) 

which spoils the numerical convergence of stochastic integration. 

The calculation of the "hard diagram" is performed with a five- dimensional nu- 
merical integration after expressing the vertex in terms of master integrals. A naive 
numerical integration of the expression in (I4.4p . without a reduction of the vertex to 
master integrals, is also five dimensional. However, our reduction method has the ad- 
vantage of removing integrable singularities which may emerge in certain kinematic 
limits. Importantly, the naive numerical integration method becomes unstable for 
n > 8, while the combination of the reduction with numerical integration is stable 
and efficient. 

As we have mentioned, the number of "master" functions required is independent 
of n. However, there are many possibilities for the kinematic invariants which enter 
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as arguments of these functions - they can be squares of either lightlike or massive 
momenta. However, we notice that evaluating the same master integrals with some 
of the kinematical invariants equal to zero may (and often does) yield a different 
structure of infrared or integrable singularities. Therefore, these cases are treated 
distinctly in our numerical approach. 

In the rest of the paper, we present explicit results for n = 6, 7, 8. This exhausts 
all possibilities for the distinct configurations in the evaluation of the diagrams of the 
previous session. The computation of n > 7 polygon Wilson loops proceeds with an 
identical algorithm as for n = 7. 



6 Six-point Wilson loops 

An analytic form of the six-point remainder function (ui, 112,113), where the six- 
point cross-ratios ^1,2,3 have been defined in (I2.25p . is currently not available. However 
we have used our numerical methods to map this function in a number of ways, as 
we will now report. 

Before doing this, we note that some numerical calculations of values of the six- 
point remainder function were presented in [23]. As an initial check on the validity 
and consistency of our methods, we compared our results with those of this refer- 
ence; we found complete agreement. Recall that we are using a different definition 
for the remainder function than [221I23], which for the hexagon Wilson loop implies 
TZ^ L (ui, U2, U3) = 7Zq hks («i, m 2 , M3) — cw- In [23] cw was found to be equal to 
12.1756 with an absolute accuracy of about 10~ 3 and we observe that this is close 
to 7r 4 /8 ~ 12.1761. For example, studying the collinear limit with our numerical 
routines we found 



in agreement with [22] . Note that one cannot simply calculate values of this function 

with one of the u variables set to zero, as the errors grow as any of the variables 
approaches zero; however one can plot the functions obtained for TZ^ L (ui,u, 1 — u) 
for various values of Ui and see that this function becomes flat as Ui — > [22J. We 
have done this with our routines and find the value ± 0.01, in agreement with 
the value 12.1756 found in [23] for TZ^ nKS in the collinear limit. Furthermore, (16.1 j) 
implies that 7^ L (0,0, 1) = 0, which is consistent with the predictions of [20 tl66l - f68] 
derived in the multi-Regge kinematics (at least in the case where all the kinematic 
invariants are defined in the Euclidean region — s ^> — — tj > 0). 

In Table 2 of [23], a number of values of 7£g HKS are also listed for different kine- 
matics. We have checked the values of the remainder function for all these inputs and 




(6.1) 
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are in perfect agreement with the quoted results up to the ubiquitous constant cw- 

The remainder function TZ$ (ui, u 2 , u 3 ) is also symmetric under permutations of 
the three cross-ratios. We have checked this in various particular cases and it is also 
amply confirmed by the results plotted in the graphs below. Before we discuss some 
plots of our numerical results in more detail we wish to add a couple of comments on 
the dual conformal invariance 

Of^ 

(ui, u 2 , U3), which, if correct, implies that in the 
six-point case the remainder function should depend on the explicit gluon momenta 
only through the cross-ratios Ui,u 2 , U3. We have confirmed this expectation by various 
numerical tests where we held the cross-ratios fixed but varied the Mandelstam vari- 
ables x 2 i+2 and x 2 i+3 . Furthermore, we have tested conformal invariance for kinematic 
points that obey the Gram determinant constraints (strictly four- dimensional kine- 
matics) and for kinematic points that do not obey the Gram determinant constraints. 
We always found perfect agreement within our numeric accuracy for TZ^^x, u 2 , 113) 
as long as the cross-ratios were held fixed. 



(ui,u 2 ,u 3 ) 




TZ^ L (B) 


ft 6 WL (C) 


(1/9,1/9,1/9) 


5.18056 


5.18096 


5.18102 


(1/4,1/4,1/4) 


1.08916 


1.08916 


1.08919 


(1,1,1) 


-2.70814 


-2.7066 


-2.70657 


(100,100,100) 


-2.09134 


-2.09204 


-2.09228 



Table 1: Checks of conformal invariance of the remainder function . In each 
horizontal line we present values of TZ^ L for different kinematic points (A), (B) and 
(C) that yield the same cross-ratios (ui,u 2 ,u 3 ) . We find that within our numerical, 
absolute errors ±0.01 the values match perfectly. Note that this estimate of the errors 
is rather conservative and that the actual error is closer to ±0.001. 



In the following we give a couple of explicit examples to demonstrate dual con- 
formal invariance of the remainder function at six points. To be more concrete we 
consider four kinematic points u\ = u 2 = u 3 = 1/9, 1/4, 1, 100, with (A) xf i+2 = —1, 

4- 



(B) x 2 3 — x 24 — 1, x 35 — x 2 m — x\ x — x\ 2 — 2 and (C) x\ 3 — x 24 — 1, x\ 



x 2 m = —2, = x\ 2 = —3. The numerical results for these kinematic points are 
collected in Table [H In general these kinematic configurations do not obey the Gram 
determinant constraint, but we have checked for numerous values of (ui,u 2 ,u 3 ) that 
TZ^ L is independent of the Gram determinant constraint as was also observed in [23] 
for one particular set of cross-ratios. We will here only discuss in detail the case 
u\ = u 2 = W3 = 100 for which one possible kinematic point, consistent with the Gram 
determinant constraint, is given by: 



2 2 2 2 20±3V42 2 2 

ry rf» T* T* '"Y* />"' I 

x 13 — x 35 — " t 46 — x 62 _ ^ ' 24 _ 51 _ ' 

x 2 -x 2 -- 1 x 2 _ 389 + 60y/42 

X\4 — x 25 — ^ , x 36 — 2Q ■ l u - z J 
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For this particular kinematic configuration the numerical evaluation of the remainder 
functions yields —2.099 with an absolute error of about ±0.02, which is in agreement 
with the last row of Table [TJ 




Figure 8: A plot of the remainder function of the hexagon Wilson loop with u\ = u<i = 
1*3 = u. For the location of the minimum of this function we found a numerical value 
ofu = 3.83 ±0.01. 

We now turn to describe additional numerical results we have found for TZ^ L (ui, u 2 , u 3 ). 
In order to explore possible analytic expressions for this function, we first considered 

F 6 (u) = nf L (u,u,u). (6.3) 

A plot of this function is given in Figure [SI Salient features are the minimum value of 
Fq{u) which is — 3.60(±0.01) at u = 3.83(±0.01) and the asymptotic value F 6 (u) — > 
-0.67(±0.05) as u — y oo. Another special value is F 6 (l) = -2.706(±0.007). It is 
interesting to observe that the minimum of Fq(u) and Fq(1) are well approximated 
by transcendentality four numbers, namely — 1= ~ —3.6077 and — |g ~ —2.7058, 
respectively. 

Some of the generic features of Fq(u) are reproduced well by transcendentality four 
functions. For example, we comment that ln(w) Lis(— l/u) + (1/3) ln(w) 2 Li2(— 1/w) 
matches the asymptotic behaviour of Fq(u), behaving as const ■ ln(-u) 2 at u — > and 
going to a constant as u — > oo, as well as matching its general shape. Whilst this 
is perhaps encouraging, at present we do not have a global match of our numerical 
results to an explicit function with transcendentality four. We also observe from our 
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data that F§{u) ~ ln 2 (w) for asymptotically small values of u. If the multiplicative 
constant is equal to rational number times 7r 2 , this is consistent with transcendentality 
four behavior. 




Figure 9: This graph contains five plots of the remainder function of the hexagon 
Wilson loop with u\ = u, ui = v and 113 = w. The cross-ratios u and v vary between 
and 10 5 while for w we have chosen five fixed values: w = 1 blue plot, w = 10 
green plot, w = 100 yellow plot, w = 1000 orange plot, and w = 10000 red plot. 

A more complete picture of the structure of the remainder function may be ob- 
tained by exhibiting a selection of slices of 7Z^ L (u, v, w) at different values of w. Due 
to the slowly changing behaviour of the function as u and v vary, it proves instructive 
to give log-based plots - in the following the u,v coordinates run from 10~ 3 to 10 5 
(listed as —3, 5 in the Figure). In Figure [9] and Figure [TOl we present from different 
viewpoints plots of the function TZ^ L (u,v,w) for these values of (u,v), in five cases 
where w = l, 10, 100, 1000, 10000. 

One may make a number of comments regarding these plots. Firstly, the symmetry 
of the function under the interchange of u and v is apparent (and is manifest in the 
actual data). Secondly, for the values of w considered, for small u, v the remainder 
function takes a large negative value for w large, and increases as w decreases. The 
order of these hyper- surfaces reverses as u or v increase, and for large values of these 
variables, the remainder function becomes increasingly negative (as of course required 
by the symmetry property and the behaviour at large u). For all three variables large, 
the remainder function approaches a constant which is equal to the asymptotic value 
of Fq(u) (about —0.67). In general, it is apparent from FigureslHlTTOlthat the remainder 
function is rather smooth for all values of the cross-ratios. 
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Figure 10: This graph contains five plots of the remainder function of the hexagon 
Wilson loop with U\ = u, u 2 = v and u 3 = w. The cross-ratios u and v vary between 
10~ 3 and 10 5 while for w we have chosen five fixed values: w = 1 blue plot, w — 10 
green plot, w = 100 yellow plot, w = 1000 orange plot, and w = 10000 red plot. 

7 Seven-point Wilson loops and collinear limits 

In this section we wish to address two separate issues. 

Firstly, we present numerical evidence that the seven-point Wilson loop remainder 
function is a function of the appropriate seven-point cross-ratios. As anticipated in 
Section 12.11 we define these cross-ratios without requiring the Gram determinant 
conditions, therefore we expect to have seven cross-ratios at seven points. As a 
basis of seven independent cross-ratios at seven points, we will choose the following 
quantities, 

«14, "25, "36, «47> «15, «26, "37 , (7-1) 

where u%j is defined in (12.231) . 

Next, we will study how the remainder function behaves under collinear limits. 
In particular, we will see that the seven cross-ratios in ( 17. ip will naturally flow into 
four parameters, three of which are naturally related to the three six-point cross- 
ratios; we will then present evidence that there is in fact no dependence on the fourth 
parameter, related to the parameter z introduced in the collinear limit (see (17. 3p 
below). As we shall discuss, our results support the conjecture that the remainder 
function of the Wilson loops should be equal to the corresponding remainder function 
on the amplitude side. 
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7.1 Seven-point remainder function and conformal invari- 
ance 

We can now compute two-loop contributions to the logarithm of the seven-point 
Wilson loop for arbitrary kinematics. There are fourteen kinematic variables formed 
by seven two-particle invariants and seven three-particle invariants. In this case we 
keep the two-particle invariants as independent inputs, and trade the three-particle 
invariants for the seven conformally invariant cross-ratios defined above. 

The sum of all relevant diagrams gives rise to the two-loop contribution to the 
logarithm of the Wilson loop. After subtracting from it the known BDS expression 
we find the remainder function . It follows from our numerical calculations that 
7£y is independent of the non-conformal input (in this case the seven two-particle 
invariants) and is only a function of the cross-ratios: 

7?wl = n^ h (u u ,u 2 5,u 3G ,u i7 ,u 15 ,u 2 e,u 37 ) . (7.2) 

This function is also invariant under cyclic permutations of all it's and under the 
reflection symmetry (which exchanges the clockwise with the anticlockwise ordering 
of w's). 

Below we give some explicit examples to demonstrate dual conformal invariance 
of the remainder function at seven points and invariance under cyclic permutations 
and reflection. To be more concrete we consider kinematic points for various values 
of the conformal cross-ratios with (A) xf i+2 = — 1 and (B) xj i+2 — — i for i — 1 . . . 7. 
The numerical results for these kinematic points are collected in Table [2J 

We have computed the remainder function for many other values of the cross- 
ratios. In Figure [11] we display the remainder function when all cross-ratios are 
equal. 



7.2 Simple collinear limits 

In this section we present numerical evidence that collinear limits of n-gon Wilson 
loops with n = 7 behave in the same way as the corresponding amplitude collinear 
limits. Collinear limits have been used very often as tools to check the consistency 
of ansatze for the expression of infinite sequences of scattering amplitudes, see for 
example [35J . In much the same way, by showing that Wilson loops, and in particular 
the remainder function, have the expected collinear limits, we can provide further 
evidence in support of the conjectured duality between the finite parts of these two 
a priori completely different quantities. 
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(Uu, U 2 5, U 3G , U47, M15, U 2 6, M37) 


K 7 (A) 


K 7 (B) 


(1,1,1,1,1,1,1) 


-3.85627 


-3.85732 


(1/4, 1/4, 1/4, 1/4, 1/4, 1/4, 1/4) 


8.13063 


8.13272 


(1/4,1,1,1/4,1,1,1) 


-4.40748 


-4.40651 


(1,1/4,1,1,1/4,1,1) 


-4.40657 


-4.40056 


(1,1,1/4,1,1,1/4,1) 


-4.40654 


-4.40559 


(1,1,1,1/4,1,1,1/4) 


-4.40746 


-4.40617 


(1,1/2,1,1,1,1/4,1) 


-4.27219 


-4.27108 


(1,1/4,1,1,1,1/2,1) 


-4.27224 


-4.27049 


(1/4,1,1/4,1,1,1,1) 


-4.63668 


-4.63696 



Table 2: Checks of conformal invariance and invariance under cyclic permutations of 
the u's and reflection symmetry of the remainder function IZ^. In each horizontal 
line we present values of1Z 7 for different kinematic points (A) and (B) that yield the 
same cross ratios. We find that within our numerical, absolute errors, which range 
between ±0.001 and ±0.01 for individual kinematic configurations, the values match 
nicely. 

We begin with a brief review of the universal behaviour of scattering amplitudes 
under simple collinear limits. In this limit, one selects two adjacent momenta p a and 
Pb, and sets 

Pa = zP , p b = (1 - z)P . (7.3) 

The collinear limit is taken by letting P 2 — > 0. Under this limit, scattering amplitudes 
behave in a well-known, universal way. Consider for instance a one-loop scattering 
amplitude, ^~ loop . When the two momenta p a andp& become collinear, the amplitude 
is known to factorise as [3514381 

(7.4) 

,{a + by,...,n) 
.,(a + b) a ,...,n) . 

Split tree are tree-level splitting amplitudes, whose explicit forms can be found, for 
instance, in [69J. Split 1- loop is a one-loop splitting amplitude. Explicit formulae for 
this one-loop splitting amplitude, valid to all orders in the dimensional regularisation 
parameter e, were presented in [70] and [71] . We quote here the result of [7T] for the 
M = 4 theory: 

Split ^ loop (a\& A& ) = Split^ c (a A %6^)r«(e^, Sab ) , (7.5) 



^- loop (l,...,a\^ 



n) > 



E 



Spiit^C^, & Ab ) ^-i op (i, - - • 

+ Split ^ loop (a\& Ab ) 
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Rj(u, u, u, u, u, u, u) 
300 r 




-300 L 
0.001 



Figure 11: A plot of the remainder function of the seven-sided Wilson loop with all 
cross-ratios equal to u. 



where, to all orders in e, 



\ c r f—s c 



v 1 ) / \ 1 



e 2 V /i 2 



F 1, -e, 1 - e. 



z-1 



F l,-e,l-e 



and cr is denned in ( 13. lip . 



2-1 
(7.6) 



We now move on to consider the behaviour of the remainder function lZ n as 
defined in ( I3.24p under collinear limits, following [21]. Consider then the two-loop 

(2) 

term J\4 n (e) in the expansion of the amplitude, and write it as 



M%\e)-\(MV{e)f = f^(e)M^(2e) + C< 2 > + TZ n + 0(e) 



(7.7) 



where f^{e) := / (2) + ff ] e + ff ] e 2 and C® are defined in (E3D, Using (17311 

and (17.51) . one sees that, under a simple collinear limit, the scalar function Ai^ must 
behave as [2] 



ofeJ 



(7.f 



(i) 



,(2), 



It was shown in pQ that splitting amplitudes obey an iterative formula identical to 
the homogeneous form of the BDS conjecture for the amplitude, i.e. 

rf{e-z,s ah )-\{r^\e-z,s ab )) 2 = f&\e)rf{2e-z, s ab ) + 0(e) . (7.9) 
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Using (I7.8P and (17. 9L one sees that, under a simple collinear limit, 

M 2) (e)-^(-Mi 1) (e)) 2 -/ (2) (e)M 1) (2e) 

" ^(-M^i(e)) 2 - /( 2 )(e)-Mi 1 2 1 (2e) . (7.10) 
Equation (17. 7p defined the finite remainder function 7£ n of the amplitude as 

72* = M%>(e) - \{M^\e)) 2 - f^(e)M^(2e) - C® + Q(e) , (7.11) 

and it follows from ( I7.10p that in the simple collinear limit TZ n — > TZ n -i, as antici- 
pated in (1X201 . 

What about simple collinear limits of Wilson loops? If the duality with amplitudes 
holds, we expect that the Wilson loop will have the same collinear limits as the 
amplitude, as discussed in (13.91) . 

Let us now specify this discussion to the seven-point Wilson loop case. Specifically, 
in the simple collinear limit of a seven-point amplitude one expects to find 

K™ h («i4, u 25 , w 3 6, «47> u 15 , u 26 , u 37 ) T$? h {u u u 2 , u 3 ) , (7-12) 

where the seven-point basis of cross-ratios is defined in (17. ip . while the six-point basis 
in (12.251) . Notice that on the right hand side of ( I7.12p we do not allow for an additional 
constant term. 

For concreteness we take p$ and p 7 to be collinear. We set 

Pe = xq — X7 = zP , p 7 = x-j — xi = (1 — z)P , (7-13) 

where as usual P 2 — > in the collinear limit. In this collinear limit, the cusp at x-j 
is therefore "flattened". This collinear limit of the seven-point kinematics is charac- 
terised by 

X 16 ~^ ' X 21 = (1 ~~ Z ) X 26 ' x 57 = ZX 51 i '-14) 

x 37 ~ V 1 z ) x 36 "+" ZX 13 i X 47 ~~ \ l z ) x m ' ZX 14 5 

with all other x?- segments unchanged. 



As one can see in Figure [T21 there are three cross-ratios which do not pass through 
the cusp at x-j. The first two are members of our basis (17. ip . whereas the third one 
is a product of members of our basis, 

2 2 2 2 2 2 

ty ry** ry^ T* ty" 

■ L \h x 3i x 26 x 35 x 13 x 46 (n -i r\ 

ry^ /y ^ /y^ /y-^ ,y-^ 

X U X 35 X 25 J '36 x W Jj 36 
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7 



7 



Figure 12: The cross-ratios Uu, U25 and u 3 qU 37 of the heptagon Wilson loop which 
do not include the cusp at x 7 . The red dashed lines depict a& segments on the right 
hand sides of (|7.15p . 

These ratios agree precisely with the three variables of the six-point case, 

(7) , (6) (7) k (6) (7) (7) k (6) ,„ a ^ 

where the superscripts on the it's denote the number of edges of the corresponding 
Wilson loop and u±,U2,u 3 are the six-point cross-ratios defined in [22] and given in 
(I2.25p . Specifically, the seven-point cross-ratios defined in (17. ip become (see (12.231) ) 

22 22 

X 24 X 15 x 26 x 35 
M14 — —9 5- , «25 — —9 q- , 

x 25 x 14 J/ 25 x 36 

•^46 ^13 I 1 ^ x 36 „, _ zx 14 „, _ n 

"36 — — 2 2 1 /1 _ — \~~2~~ ' " 47 ~~ — 2 _i_ /1 _ — TIT ' 15 ~~ U ' 

2*3g %%u t (1 z ) x ^q 14 ~r z )^4& 

2 ^,2 ,^2 1 /1 



v l ^)^36 _ X 13 2X 14 + (1 2^)^46 ---1 -1 

^26— 2 T /1 _ — \~~2~~ ' M 37 — — 2 2 , /1 _ — ■ K'- 1 ') 

zx 13 + I 1 Z ) X 3G X 14 zx 13 + I 1 Z ) X 36 



As expected, it follows immediately that uu, U25 and u 36 u 37 are equal to the three 
cross-ratios of the six-point case, as dictated by (17. 16ft . In addition, (17. 17ft imply that 
there are three constraints on the remaining four variables, 

«15 = , U A7 + U 26 M 36 = 1 , «26 + ^37 M47 = 1 • (7-18) 

Taking this into account, and solving the two constraints in (I7.18P for m 47 and u 2 6 
we conclude that the collinear limit relates the seven-point remainder function to the 
six-point one as follows: 

T>WL/ \ -nWL ( 1—^36 1 — U 37 
TZ 7 {Ui4,U25,U 36 ,U47,Ui 5 ,U26,U 37 ) ->■ K 7 [ Uu, U 2 5, «36, Z , 0, , U 3 7 



1 - u 37 u 36 1 - u 37 u 36 
T^ L (u 37 u 36 ,uu,U2 5 ) ■ (7.19) 
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Note that three of the seven variables on the left hand side of ( 17. 191) are constrained, 
leaving four variables free. Since the right hand side is a function of only three vari- 
ables, it follows that the left hand side actually does not depend on one combination: 
M37/M36 := k. Thus we can rewrite ( I7.19P in terms of the six-point cross-ratios Ui and 

Kj cLS 



1iy h u 2 , u 3 , VW^, 1 _ » °> 1 _ 1 » v 7 ^ = u 2 , u 3 ) , (7.20) 

and note that the left hand side must therefore be independent of the variable k. 
This can be thought of as the z-independence of TZ^ h in the collinear limit - which 
is precisely the feature one would expect from the scattering amplitude. 

We have computed TZ^ L in the collinear limit and have confirmed that (I7.20p 
does hold within the errors and that no dependence on k is found (see Table E] and 
Figure [T31 for example). Again, we stress the absence on the right hand side of (I7.20p 
of any additional constant term. 



{u u u 2 , Uz) 


ti^{k = 2.5) 


7^(« = 4.9) 


-nWL 


(1/10,1,1) 


-2.78972 


-2.76053 


-2.73441 



Table 3: Checks of the collinear limit TZf L ->■ TZf L . We present U^(k) := 

(u2, u 3 , a/ Ui/ k, 1 0.01, 1 y/uiKj for different values of K, together 

with its collinear limit TZ 6 (ui, u 2 , u 3 ) for (ui,u 2 ,u 3 ) = (1/10, 1, 1). Within numerical 
errors the values agree. 

In the following section we compare these results with what can be learned from 
the multi-collinear limits. 



7.3 Multi-collinear limits 



Here we would like to derive the multi-collinear equivalent of the general reduction 
formulae in (13. 9p . (17.19p - (17.20p for n-gon Wilson loops. 

The first non-trivial case is a triple collinear limit of a six-point configuration 
considered in Section 5 of jUJ. In the limit where p±, p$ and p§ become collinear one 
has 

p 4 : = X4-X5 = z 1 P , p 5 = x 5 -x e = z 2 P , p 6 = x 6 -x 1 = z 3 P , z 1 + z 2 + z 3 = l. 

(7.21) 
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nWL „WL 



-0.2 
-0.4 

-0.6 
-0.8 
-1.0 
-1.2 
-1.4 




1.0 



Figure 13: A plot of the seven-point remainder function minus the six-point remainder 
function in the collinear limit. The kinematics are given in (I7.20p with (ui,u 2 ,u 3 ) = 
(1/10, 1, 1) but with a non-vanishing ui§ which varies along the x-axis (ui§ vanishes 
in the collinear limit). For the blue graph we have k = 0.25 and for the purple 
graph n = 0.49. We see that in both cases the difference tends to zero with m 15; 
confirming (I7.20p . 



In this limit, the two-loop scalar function Ai^ behaves as [21] 

\S456 S456 / V S456 S456 



(7.22) 



where and are the one- and two-loop triple splitting amplitudes. The two- 
loop triple splitting amplitude does not satisfy an iteration relation similar to the 
that found in [TJ for the simple splitting amplitude (17. 9p |21j . It can however be 
decomposed into a term with does, plus a two-loop finite remainder, as 



,(2) 



r 



(2)BDS 

s 



+ K . 



(7.23) 



It was shown in [21] that 1Z is nothing but the finite remainder function at six points 
evaluated in the triple collinear kinematics. Indeed, in the triple collinear limit the 
six-point cross-ratios remain independent and do not vanish, — > Ui, i = 1,2,3, 
where 



1 



Ui 



S45 



1 



«2 



S56 



«3 



ZlZ 3 



Hence, 



1 — Z3 S456 1 — Z\ S456 (1 — Z\){\ — Z 3 ) 

1Z 6 (ui,U2,u 3 ) -> TZe(ui, u 2 , u 3 ) . 



(7.24) 
(7.25) 



36 



If we now take the triple collinear limit of the remainder function as defined from 
(17TT|) . and use (172211 and (172311 . we get 

7^ 6 (mi,M2,%) -> Tli + 7^(wi,«2,«3) • (7-26) 
Using that 7^ 4 = 0, and comparing (17.261) to (I7.25P we get at once that 

TZ(u 1 ,u 2 ,u 3 ) = 1Z 6 (ux,U2,u 3 ) . (7.27) 

Assuming further dual conformal invariance, the conclusion of [21] is that determining 
the remainder function TZ e in the triple collinear limit is equivalent to determining it 
in full generality. 

We can now derive the triple collinear limit of the two-loop n-point remainder 
function for general n > 6. Using (I7.27P we have 

TZ n -> K n _ 2 + 7£ 6 (ui,U2,w 3 ) , (7.28) 

where (ui, «2> u 3 ) are still defined by (17.241) (in the case where the collinear momenta 
are p 4 , p 5 and p 6 ). 

What about Wilson loops? Similarly to our discussion of simple collinear limits, we 
expect that, for Wilson loops, the triple collinear limit of the Wilson loop remainder 
function is given by 

TZ^ L (Mi4,M25,M 3 6,M47,Mi5,M 26 ,M 37 ) -> TZ^} 2 + TZ^ L (lij , u 2 , u 3 ) . (7.29) 

Now consider a triple collinear limit of a heptagon Wilson loop, where p$, p$ and p-j 
are collinear, 

p 5 = x 5 -x 6 = z x P , p 6 = x 6 -x 7 = z 2 P , p 7 = x 7 —xi = z 3 P , z 1 + z 2 + z 3 = l. 

(7.30) 

In this limit, 

^wl ^ n™ L (u u u 2 ,u 3 ) , (7.31) 
since there is no five-point remainder function. 

For our present case of collinear p 5 , p 6 and p-j the variables Ui read 

l — ^s 5 2 1-Z X X\ % ' 3 (1 -«i)(l -Z 3 ) ' 

The triple collinear limit ( I7.30P of the seven-point kinematics gives 

£?5 ~ A& ~ ^57 -> > 
2 2 2 /i \ 2 2 2 2 \ 2 

^36 = (1 - *l)&35 + Z l X ls , ^37 = ( X ~ Z s) x l 3 + z 3 x l 5 , ( 7 - 33 ) 
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with all other x\^ segments unmodified. The cross-ratios Uij take the form: 

(1 " z l) x 35 
U\4 — U , U 25 — 2 T~7i \~~2~ 

ZxXlz + (1 - z 1 )x^ b 

2 1 rr2 



Z l Z 3 X 35 + (! - ^)a;?3 1 



,r 



57 „, x ^16 



1 - 2 3 ZiXf 3 + (1 - 1 - 2 3 xf 5 1 - Z x xf 5 

^3 + (1 - Zx)x 2 35 (1 -Z 3 )X 2 13 



U2& ~ M37 = — (7 34) 

These relations imply two things. First we note that there are three constraints on 
the seven variables, 

u u = , u 25 + u 36 u 37 = 1 , u 31 + u 25 u 26 = 1 , (7.35) 

leaving four variables unconstrained, out of which the three conformal ratios coincide 
with the Ui variables of (17.321) . 

M47 = ui , M15 = u 2 , u 26 u 36 = u 3 . (7.36) 

Taking this into account, and solving the two constraints in (I7.35P for Uu and M35 we 
conclude that the triple collinear limit relates the seven-point remainder function to 
the six-point one as follows: 

»WL/ n 1 _ n 36 „, „, l-«26 \ _ ^\\L, 



Tl 7 0,- , «36,«47, «15, «26, " =^6 ( M 47, «15, ^26^36) • (7.37) 

\ l-u 2Q u 3G l-u 36 u 26 J 

Notice that since TZ^ L is invariant under cyclic interchange of its variables, this is 
exactly the same equation as we had before in the simple collinear limit ( I7.19P (alas 
expressed in terms of ii 36 , M47, u i5 and u 26 ). 

We conclude that the simple and the triple collinear limits described by (I7.19P and 
(I7.37P give identical information about TZ^ L . As before, the left hand side of (I7.37P 
cannot depend on one particular combination of cross-ratios, namely u 3 q/u 2 q := k. 
Thus, we can rewrite (I7.37P as 

,\vi, ; , 1 t=— - - f—JZ 1 - \ WL 



TZ 7 0, — = — ,^/u 3 ~k,u 1 ,u 2 ,\/u 3 /k, — = =^ 6 (m,u 2 ,u 3 ) , (7.38) 

y l-u 3 1 - u 3 J 

and note that the left hand side must be independent of the variable k. 

Finally, we have investigated the quadruple collinear limit, which is the highest 
non-trivial multi-collinear limit one can take on the seven-point kinematics. We have 
found in this limit that all seven cross-ratios are (a) mutually independent, and (b) are 
expressed entirely in terms of the multi-collinear kinematics (i.e. they are functions 
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of zi, . . . , Z4 and ratios of kinematic invariants involving only the collinear momenta). 
In this way, the quadruple collinear limit does not add any non-trivial functional 
constraints on the TZ-/, however, it elucidates its physical (scattering amplitudes- 
based) meaning, 

TZf L (u lAj U25,u m ,u 47j u l5 ,U2 6j u 37 ) = Asplit 4 := rgj - rfl BDS , (7.39) 

where Asplit 4 is the normalised two-loop level part of the quadruple splitting function 
r^l which is not already accounted by the BDS contribution r^ 2 4 BDS . Hence, similarly 
to the six-point case discussed in previously [2TJ, we see that the remainder function 
7^-7 is entirely determined by the quadruple splitting function. 



8 Eight-point Wilson loops and beyond 



It is natural to seek beyond the encouraging results at seven points given above, and 
see if these persist at eight points. The discussion earlier indicates that there are 
twelve independent conformal cross-ratios in this case (we take the external momenta 
to be on shell and do not impose the Gram determinant constraint), whereas there 
are twenty independent momentum invariants. These independent invariants may be 
taken to be 

2 2 2 2 2 • 1 A (O 1 \ 

X i i+2i X i+i i+6' X i i+3> X i+i i+7i X ii+A' 2 — 1, . . . , 4 . \P- L ) 

We will use the following twelve cross-ratios: 

u ii+3 , i = l,...,8, u ii+ 4 , z = l, ...,4, (8.2) 

and label these Ui, ...,ui2- Instead of the twenty momentum invariants given in (18. ip . 
we will use these twelve cross-ratios, plus the following eight momentum invariants 

X i+5 i+8-i X ii+i i ^ = 1, • • • , 4 , (8-3) 

which we will call mi, . . . The remaining momentum invariants xf i+2 , xf +4 i+6 

and xf +1 j +4 , for i — 1, . . . , 4, are then dependent variables. 

The first question to study is whether the eight-point Wilson loop remainder 
function TZ^ L is only a function of the twelve conformal cross-ratios (18. 2p . and not 
of the additional eight invariants (18.31) . To do this, one may fix a choice of the 
cross-ratios, then calculate the eight-point remainder function for various choices of 
kinematics, corresponding to different choices of the variables (18. 3p . For example, in 
Table H] we have listed some numerical results for the case of all the cross-ratios of 
(18. 2p equal to one. We find similar results for more generic values of the cross-ratios 
- an example is given in Table 
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(mi, . . . ,m 8 ) 


ft 8 WL 


(-1,-1,-1,-1,-1,-1,-1,-1) 


-4.603 


(-2,-2,-2,-2,-2,-2,-2,-2) 


-4.602 


(-1,-2,-4,-8,-1,-2,-4,-8) 


-4.605 


(-5,-3,-5,-3,-1,-3,-5,-7) 


-4.605 



Table 4: The remainder function 1Z% for Ui — — • ■ • — U\i = 1 and different 
choices of the other independent invariants ( 18.31) . The errors in lZ™ h are approxi- 
mately 0.02. 



(mi, . . . ,m 8 ) 




(-2,-3,-4,-1,-5,-6,-7,-8) 


5.993 


(-1/3, -1/4, -1/9, -1/2, -1/8, -1/7, -1/6, -1) 


5.984 



Table 5: The remainder function TZ^ L for the choice of cross-ratios («i, . . . , U12) = 
(2,3,4,1/2,1/3,1/4,1/5,1,1/5,1/6,1/7,1/8) and two choices of the other indepen- 
dent invariants (18. 3p . The errors in 1Z$ are approximately 0.04. 

These results support the conjecture that the eight-point remainder function 7£g 
is a function of the twelve cross-ratios ( 18.21) alone, and not of any other independent 
additional momentum invariants. 

Given this, a second question concerns symmetries of the function lZ^ L (ui, Uu). 
The Wilson loop is invariant under cyclic permutations of the external momenta as 
well as under parity. For the case of eight points, this implies that the remainder func- 
tion TZ^ h (ui, W12) should be invariant under cyclic permutations of the first eight 
and last four cross-ratios simultaneously, as well as being invariant under the simul- 
taneous reversal of both the first eight and last four cross-ratios. We find numerical 
agreement with this - for example, 

^ L (2, 2, 1,1, 1,1, 1,1, 1,1, 1,1) = -3.712, 

TZ^ L {1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1) = -3.712, (8.4) 

with errors ~ 0.02. We have also seen numerically in a number of cases that there is 
no invariance under more general permutations of the cross-ratios. 

Let us now consider collinear limits of the eight-point remainder function. As an 
example we take pi and p$ to be collinear. We set 

p 7 = x 7 - x$ = zP , p 8 = x s - Xx = (1 - z)P , (8.5) 

where as usual P 2 — > in the collinear limit. In this collinear limit, the cusp at xg is 
"flattened" . In this limit we find that the eight-point cross-ratios (uu, U25, M36, M15, «26) 
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reduce directly to the cross-ratios in the seven-point case with the same names, 
Ui6 — > and the seven-point cross-ratios (1*47, W37) are given in terms of the eight- 
point ones as {u4 7 U4 8 ,u 37 u 38 ) respectively. Finally, we have the following relations 
amongst the eight-point cross-ratios: 

U27U37U38 = -1 + U 2 7 + u 38 , 

(-1 + U 27 + U 38 )U A7 = U 38 (l - M 58 ), 
^48^58 = 1 - ^27^37, 

U 37 (~U 38 - U47 + U 38 U A7 + U 38 U 58 ) = -1 + «58, (8-6) 

which are solved by 

-1 + u 3S 



u 2 7 



U58 



-1 + U37-U38 

-1 + U 37 U 38 + ^37^47 ~ ^37^38^47 
-1 + U 37 U 38 



M 4 8 = —— • (8.7) 

-1 + M37M38 + M37M47 - U 37 U 38 U A7 

The three variables (u 37 ,u 38 ,U4 7 ) in the above are then freely specifiable. 

Analysis of the remainder function given earlier implies that in the collinear limit 

Hence in the collinear limit (18. 5p one should have 

^("14. M 25, «36, M 47, «58, «16, U 27 , U 38 , U 15 , U 26 , U 37 , U 48 ) -» 
TZ^ (U U , U 25 , «36, M 4 7, U* 58 , 0, U* 27 , U 38 , U15, «26, U 37 , U* 48 ) 

= 7?.^ L (Mi4, «25, «36, U47U*4 8 , «15, «26, ^37^38) , (8.9) 

where the stars in the above indicate that the solutions (18 .7\i are to be inserted. 

One can test this directly. For example, for the choices of values for the indepen- 
dently specifiable eight-point cross-ratio variables 

(U U , «25, «36, "47, u 83 , u 15 , U 26 , u 37 ) = (1, 1, 1, 1, 1/2, 1, 1, 1/2), (8.10) 

and taking w 6 i = 0.001 one finds 

7^ L (wi4, U 2 5, «36, U47, u* 58 , 0.001, u* 72} U 83 , U 15 , M 2 6, «37, u* 4s ) = -4.2756, 
TZf L (u u , U25, u 36 , u A7 u A8n u 15 , u 2 6, u 37 u 38 ) = -4.2906 , (8.11) 
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with errors of 0.147 and 0.005 respectively. 

The success of the above numerical tests of the conformal symmetry, functional 
symmetries and collinear limits of the eight-point remainder function supports the 
conjecture that the Wilson loop is correctly reproducing the physical amplitude at 
this level. 

There are no further conceptual or computational obstacles to generalising the 
above work beyond eight-point Wilson loops, apart from the question of the computer 
time required to numerically calculate the integrals - we stress that no new integrals 
arise in the Wilson loop approach to n-point two-loop diagrams for any n, apart from 
those which we have already discussed, and we are able to calculate all the diagrams 
introduced for generic values of the momentum variables Qi (numerically, and in a 
number of cases, analytically). 

This means that we have full numerical control over two-loop n-gon Wilson loops 
and, if the correspondence with amplitudes continues to hold, over n-point MHV 
amplitudes at arbitrary n in the planar Af = 4 theory. This should be contrasted 
with the situation where one calculates amplitudes directly. 
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A A note on conventions 



Wilson loops are computed in dimensional reduction in D = 4 — 2e uv dimensions 
with e uv > to regularise the UV divergences. To facilitate the comparison with 
scattering amplitudes (which require infrared regularisation) we introduce 

e = -e uv . (A.l) 

The perturbative expansion of the Wilson loop is characterised by (13. 2p . (13. 3ft : 

oo oo 

(W[C n }) = l + J2a l W® =«pj> l u;g>, (A .2) 
i=i i=i 

= W^-\{W^)\ (A.3) 
and (I3.2ip defines the Laurent expansion in e for the two-loop contribution, 

q \ -26 ,„„(2) (2). 



n / 2 \ e W v ' 

£ ("^) (^ + ! T i ) + ^ + °M • < A - 4 > 



(2) 

where Fn is the finite part of the Wilson loop. 

In [23] the exponent of the Wilson loop in II A. 21) was defined as aw^+2a 2 w^ + - ■ ■ , 
and thus a corresponding factor of 1/2 would need to be introduced in front of the right 
hand side in (IA.3j) and in (1A.4|) if we were to switch to their conventions. Therefore 
our singular terms w_2, w~i and the finite part Fn in ( 1A.4|) are related to the A_ 2 , 
A_i and A contributions of [25] as follows: 

„_ 2 = £ ^ , «,_! = - J] ^ , = 2 E ^ • (A.5) 

a a a 

The minus sign in the second equation is due to flA.li . 



B Most general hard diagram 



Below is the result for the diagram where a three-point vertex is attached to three 
lightlike momenta pi, p 2 , ps of the Wilson loop, which we call the "hard diagram" as it 
is the most difficult to evaluate analytically in general. These momenta are separated 
by the three, not necessarily lightlike momenta, Q 3 , Qi, Q 2 , where Q 3 is between p x 
and P2 and so on (see Figure |2]). Momentum conservation is then X/i=i(P« + Qi) = 0- 
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We also set D = 4 — 2e uv 
is considered later. 



4 + 2e where e t 



-e > 0. The special four-point case 



We write this diagram in the most general configuration a^| 

fH(pi,P2,P3;Qi,Q2,Qs) 



r l 



on) \a\_a 2 a 3 ) 



" e UV . 



1=1 



1=1 



i=l 



V 2- 



2c, 



where 



and 



-aia 2 {z x - z 2 ) 2 - a 2 a 3 (z 2 - z 3 ) 2 - aia 3 (zi - z 3 ) 2 , 



(z x -z 2 ) 2 = Ql + 2(p 1 p 2 )(l-T 1 )T 2 + 2(Q 3 p 1 )(l-T 1 ) + 2(Q 3 p 2 )r 2 
(z 2 -z 3 ) 2 = Ql + 2{p 2 p 3 ){l-T 2 )T 3 + 2{Q l p 2 ){l-T 2 ) + 2{Q l p 3 )T 3 
(z 3 - Zl ) 2 = Ql + 2(p 3 p 1 ){l-T 3 )r 1 + 2(Q 2 p 3 )(l-r 3 ) + 2(Q 2 p 1 )r 1 

The original expressions for the Z{ — Zi + \ are 

Zi-z i+1 = Q i+2 + pi(l - n) + p l+1 r i+1 , i = l,2,3. 



(B.l) 
(B.2) 

(B.3) 



(B.4) 



The expression for the numerator M has two kinds of terms. The first three lines 
involve r and a parameters, whereas the remaining three lines involve only the r 
parameters. It is given by 

M 



= 2(p 1 p 2 )(p 1 p 3 ) 


ai« 2 (l 


- ri) + a 3 axT X 




+ 2(p 1 p 2 )(p 2 p 3 ) 


a 2 a 3 (l 


- t 2 ) + aia 2 r 2 




+ 2(p 1 p 3 )(p 2 p 3 ) 


a 3 «i(l 


- r 3 ) + a 2 a 3 r 3 




+ 2«i« 2 


2{piP 2 )(p 3 Q3) 


- (PaP3)(PlQs) 




+ 2a 2 a 3 


2{p 2 p 3 ){piQi) 


- (PaPi)(PaQi) 


- (PiP2)(PsQi) 


+ 2a 3 ai 


2(p 3 Pl)(P2Q2) 


- (PlP2)(P3Q2) 


- (P2P3)(piQ2) 



(B.5) 



B.l Four-point case 

The four-point case can be obtained by setting 

Qs = Qi = 0, Q 2 = p 4 = -( Pl + p 2 + p 3 ) , (B.6) 

9 We remind the reader that we will always suppress the common prefactor defined in (14. ip from 
the expression of all diagrams. 
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where now Q\ 



0. The expression for M in ( IB. 51) then simplifies to 

M = 2(pip 2 )(piZ>3) (1 - Ti)ai(a 2 - a 3 ) 

+ 2(p 1 p 2 )(p 2 P3) a 2 a 3 (l - r 2 ) + cti(2a 3 + a 2 r 2 ) 
+ 2(p 1 p 3 )(p 2 p 3 ) r 3 a 3 (a 2 - ai) . 

In this special case we have 

(zi - z 2 ) 2 = 2(pip 2 )(l - ri)r 2 , 

(z 2 -z 3 ) 2 = 2(p 2 p 3 )(l - r 2 )r 3 , 

(* 3 - ^i) 2 = 2(papi)(l - r 3 )n + 2(pap 4 )(l - r 3 ) + 2(pu? 4 )r 1 , 
where we can set 



s 
£ 
it 



2(piP 2 ) 
2(p 2 Ps) 
2(piP3) 



2(p 3 P4) 
2faiP 4 ) 

2(P2P 4 ) 



(B.7) 
(B.8) 



(B.9) 



and s + £ 



it 



0. 



The denominator in the four-point case then simplifies to 
V := saia 2 r 2 (l - n) + £a 2 a 3 r 3 (l - r 2 ) + ct 3 ai s(l - ri)(l - r 3 ) + £tiT 3 . (B.10) 



Notice that in the six-point case there is a new diagram where all the Q^s are 
made of a single lightlike momentum; in this diagram Q\ = Q\ = Q\ = 0. One should 
find that the result for the corresponding integral X is finite in four dimensions. 



C Curtain diagram 



We call "curtain" diagrams those diagrams where two propagators connect three 
different edges, as depicted in Figure [31 We use the same notation as in the three- 
point vertex case above, with propagators stretching from p\ to p 2 (with end-points 
Z\{p\) and z 2 {t 2 )) and from p 1 to p 3 (with end-points Zi(ti) and z 3 (r 3 )). 



Just from looking at the diagram we should have the following symmetry: p 2 -H- p 3 
Q2 «-> Qz- We have 



z i{ T i) = pin 

zi(ci) = Pi o-i 

32 fa) = P2^2 +Pl + Q3 

33 fa) = -P37-3 - Q2 ■ 



(C.l) 

(C.2) 
(C.3) 
(C.4) 
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The exponentiation theorem says we should only consider the diagram where the 
internal gluon propagators cross - this gives the constraint T\ > a%. 



The diagram represents the following contribution to the Wilson loop: 



(P1P2) (pm) 



( - [z x {*,) - z 2 {r 2 )Y) 1 -^ ( - [ Zl {n) - zsfo)]*) 1 -*™ ' 

(C.5) 

Putting in the values for the end-points we get the following integral representation: 



2 Jo Vjx 7 J [ _ 2 (p 1 Q 3 )(l - - 2(p 1 p 2 )(l - tn)^ - 2(p 2 Q 3 )r 2 - Q§] 

x (P1P3) 

[ - 2(p 1 Q 2 )r 1 - 2( Pi p 3 )t 1 t 3 - 2(p 3 Q 2 )r 3 - Qf] 1_euv ' 

(C.6) 

A more symmetrical way to write this is to send a\ — > 1 — <j\ in which case the 
constraint t\ > <J\ becomes T\ + <J\ > 1 and the integrand would be manifestly 
symmetric under p 2 -H- ps, Q 2 -H- Q3, t\ -H- 01 and r 2 -H- r 3 . We have performed this 
change of variables to obtain ( 14. 6p . 



D Cross diagram (involving two sides) 



This diagram consists of two gluon propagators, stretching from sides p\ to p 2 with 
sides Qi and Q2 between, as represented in Figure HI 

The end-points of the first propagator are z\{t\) and z 2 (t 2 ), and of the second are 
Zi(ai) and z 2 (a 2 ), with 

ZiM = V\°\ z 2 (a 2 ) = -p 2 a 2 - Q 2 

Zi(ri) = Pin z 2 (t 2 ) = -p 2 r 2 - Q 2 . (D.l) 

In order to ensure the crossing of the propagators we require T\ < o\ and r 2 > a 2 . 
The diagram then represents the integral 

— / daidTo [ dri f dao (PiPz) (P1P2) 

2 Jo Jo Jo ( - fofa) - z 2 {a 2 )Y) 1 -™ ( - fain) - z 2 {t 2 )Y) 1 -™ 

(D.2) 



l-et 
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Putting in the values of the end-points gives the integral 

^ da\dr 2 [ dr\ [ da 2 ~' 



2 ./ Jo Jo f - 9fT),n„W^„ - ^T)^^, - 9n„D^„ - fl^ 1 £OT 



( - 2(j»ip 2 )cr 1 (X2 - 2p 1 Q 2 a 1 - 2p 2 Q 2 (J 2 ~ Ql) 

(pm) 

( - 2(pip 2 )TiT 2 - 2p 1 Q 2 r 1 - 2p2Q 2 T 2 - Qf) 

(D.3) 



E Y diagram 



This diagram consists of three gluon propagators, meeting at a vertex, with two 
propagators ending on side p\ and the third ending on p 2 with sides Qi and Q2 
between. The contribution of this diagram is 



Pi -P2 '' 1 



dT\dr-< 







2G(z 1 (t 1 ),z 1 (t 1 ),z 2 (t 2 )) -G(z 1 (0),z 1 (t 1 ),z 2 (t 2 )) -G(z 1 (1),z 1 (t 1 ),z 2 (t 2 )) 

(E.l) 



where z\{ti) = piTi is a point ending on the edge p\ and z 2 (r 2 ) = —Qi — p 2 T~2 a point 
on edge p 2 and 

T(l - 2e uv ) f 1 Jn Jn ^ (aia 2 a 3 )" euv 5(l - a x - a 2 - a 3 ) 



G(z l ,z 2 ,z 3 ) = ^j- — ^Hll f da x da 2 da 
T 2 1 - e uv Jo 



-aia 2 zi2 - aia 3 zi3 - a2a 3 z 23 ) 1 2euv 

(E.2) 

1 m-2e uv ) [ l a -^(l- a y^ 

d(T ^—J2 71 _N_2 M-2, T „, 



e uv P(l - euv ) 7 - (1 - a)^) 1 " 2 ^ 

where the final equality is valid whenever zf 2 = (as is the case here) and can be 
shown by changing variables to p, a with ax = (1 — p)cr, a 2 — (1 — p)(l — a), a 3 = p. 

Now it turns out that the first term in (IE. 1[) is precisely canceled by half of the 
self-energy correction to the propagator between sides p\ and p 2 with the other half 
canceling the upside-down Y diagram. We also notice that it is the combination 
of these two contributions which has the expected maximal transcendentality. The 
explicit expression for the self-energy correction to the gluon propagator in M = 4 
SYM can be found, for example, in [72]. 
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Thus, neglecting the first term, the result of the Y diagram integral is 

Pi-p 2 1 r(l-2e uv ) f 1 f l / ^ £uv (l-^)" euv 

8 e uv P(l - e uv ) J a J n T2 \ (-a(Q 1 +p 2 T 2 y-(l-a)(Q 1 +p 1 T 1 +p 2 T 2 yy- 2 ^ 



(1-a)- 



r-«UVn _ rr - £ UV 



)l-2e uv ) 



(-<j(- P2 (i - r 2 ) - Q 2 y - (i - ^(Qi + pit! + P2 r 2 yy 

Pl -p 2 1 r(l-2e uv ) [\ f\ J / ^uv^.^-euv 



y rfo" y dridr 2 (^ — 



8 e uv r 2 (l-e uv ) 7 7 Z V (_(l_ (7 )(Q 1+ p 2 r 2 )2- (7 (Q 1 +j9 1 r 1 +p 2 r 2 ) 2 ) 1 - 2 ^v 

(7 _euv (l - (r) _euv \ 

_ (-(1 - <r)(p 2 T 2 + Q2Y ~ (?(Q2 + PlTl + p 2 T 2 ) 2 y-^W ) 

_ Pl ■ P2 1 r(l - 2 £uv ) [i f 1 ^ , a~^(l - a)~^ 



da J dr\dr 2 ^ — 



e uv P(l - e uv ) y J V (-Q? - 2(Q iP2 )t 2 - <rr 1 (2(Q 1 p 1 ) + 2(p 1 p 2 )r 2 ))i-*uv 



(-(Q| + 2(Q 2 p 2 )r 2 + crr 1 (2(Q 2 p 1 ) + 2(p lP2 )r 2 )) 



1=2^) ' ( E - 4 ) 



where to obtain the second equality we have used the change of variables a — > 1 — a, 
Ti — > 1 — Tj. The final answer is manifestly symmetric under Qi -H- Q 2 . 

All two-loop diagrams are given by the above integrals, for various values of the 
momenta p iy Qi. 
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